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NONLINEAR EVOLUTION OF A DISTURBANCE IN AN UNBOUNDED 
VISCOUS FLUID WITH UNIFORM SHEAR 
by Robert G. Deissler and Burt M. Rosenbaum 
Lewis Research Center 

SUMMARY 

The evolution of a disturbance in the presence of a uniform mean velocity gradient 
is calculated by a power-series solution of the incompressible Navier-Stokes equations. 
Terms through those in time cubed are retained in the solution. For the initial condition 
a three-dimensional cosine distribution with two harmonic terms is assumed. The non- 
linear interaction of these harmonic terms produces new harmonics which in turn inter- 
act. For large velocity gradients the energy of the disturbance can grow with time. 

The results may shed some light on the maintenance or growth of turbulence in a shear 
flow. 


INTRODUCTION 

The history of a disturbance in a viscous fluid has been calculated to various degrees 
of approximation for several initial conditions in references 1 to 3. The results showed 
that the evolution of the disturbance is modified by the generation of new harmonics 
through the interaction of those already present. The effect of chemical reaction on a 
fluid disturbance was studied in reference 4. Various numerical aspects of Fourier- 
type solutions of the Navier-Stokes initial value problem were considered in reference 5. 

In the present report the results of reference 3 are extended to include the effect of 
a uniform velocity gradient on the evolution of a fluid disturbance. Thus, in addition to 
mode-mode interactions (interactions between the Fourier components of the wave), in- 
teractions between the velocity fluctuations and a mean gradient are considered. As in 
reference 3 the initial condition for the velocity disturbance is taken to consist of two 
waves with different wave number and intensity vectors. It is noted that turbulence can 
be regarded as being made up of a very large number of disturbances of the type con- 
sidered in this report, so that the results may shed some light on turbulent shear flow, 


particularly in- the early stages of development. In particular, it would seem that non- 
linear effects might be more easily investigated in the present problem than in turbu- 
lence. As will be mentioned in the section RESULTS AND DISCUSSION, there are, how- 
ever, important differences between the results obtained for simple regular disturbances 
and those for a fully developed random turbulence. The case of turbulence with a uni- 
form shear, but without nonlinear effects, is studied in reference 6. 

For simplicity a Taylor power series in time rather than the exponential method of 
reference 3 is used. The basic equations and their analysis are considered in the next 
section. 


ANALYSIS 

The following problem is considered: given the initial velocity distribution in an 
unbounded viscous fluid with a uniform mean velocity gradient, predict the motion at 
later times. For a viscous fluid with constant properties, the three-dimensional equa- 
tions of motion in dimensionless form are 
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where equation (2) is obtained by taking the divergence of equation (1) and applying the 
continuity equation, and 
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The quantity u* is a velocity component, x| is a space coordinate, Xq is a character- 
istic length, t* is the time, p is the constant density, and v is the constant kinematic 
viscosity. Note that the stars on dimensional quantities are omitted for corresponding 
dimensionless quantities. The subscripts in equations (1) and (2) can take on the values 
1, 2, and 3, and a repeated subscript in a term indicates a summation. (Symbols are 
defined in the appendix. ) 

One can break the dimensionless instantaneous velocities into mean and spatially 
fluctuating components; thus, set u. = U- + u^. Equations (1) and (2) become 
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If the mean velocity is assumed to be in the x^-direction, the left side of equa- 
tion (3) becomes [0/0t + U^(0/0x^)]u^. In this analysis we let the observer move with 
the mean velocity at each point, since we are interested in the changes with time of fluc- 
tuation levels from that point of view. Also, the mean velocity gradient is taken to be in 
the Xg-direction, constant in time and space, and equal to dUj/dxg = S. Coordinates 
Xj relative to the moving observer are given by the transformation x. = x^ - ^jUjCx^t. 
Then 
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Also, 
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Equations (3) and (4) become, in the new coordinate system, 
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where 6. - is the Kronecker delta. In addition, the equation of continuity which had the 
form 
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For the initial condition on the velocity fluctuation at time t = 0, take 
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where m' can, in the general case, be any positive integer, and a^ and q- are, 
respectively, intensity and wave-number vectors. In the actual numerical calculations 
m' will be set equal to 1, 2. Two is, of course, the smallest number of terms that can 
be retained if there are to be mode-mode interactions. 
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To get the evolution of the velocity fluctuation, a Taylor series in time- is used and 
the required initial time derivatives are obtained from equations (5), (6), and (8). The 
Taylor series is 



The quantity ^u.j in this equation is obtained from equation (8), which can be more 
conveniently written in complex notation as 

m 

where m takes on both positive and negative values (i. e. , m = ±m'), and qr m = -q! 11 
and a: m = a™. In the case considered in this report, where m' = 1, 2, the values m 
takes on are given by m = -2, -1, 1, 2. Also, the equation of continuity yields the rela- 
tion 


a^qf 1 = 0 (11) 

(Note that, although a repeated subscript indicates a summation, summations over 
superscripts are only to be carried out when explicitly indicated by the summation sign. ) 

To obtain ^3u./3tj in equation (9), equation (10) is first substituted into equa- 
tion (6) evaluated at the initial time. This gives 
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where m and n are independently assigned the values assigned to m. Inspection of 
equation (12) shows that Pq has the form 


P 0 = 



m iq m - x 
it e n 





m +q n )* x 


(13) 


5 



Substituting equations (13) into equation (12) and equating coefficients of like powers 
of e give 


0 m„m 
I™ = S 2 1 


(14a) 


m „nn m 

» mn = -I L (14b) 

4 L m , „n\/ m , n\ 

(?) q j A q i + q j) 

where equation (11) has been used to simplify the expression for Jr™ 11 . Substituting 
equations (10) and (13) into equation (5) evaluated at t = 0 gives for the first time deriv- 
ative in equation (9) 
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and a? 111 has been written so as to maintain symmetry in the superscripts m and n. 

To obtain ^u^/st^ in equation (9), differentiate equations (5) and (6) with respect 

to time. This gives at t = 0 
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where (9p/3t)Q is given by 
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Substituting equations (10), (13), and (15) into equations (18) and (17) and proceeding as 
was done to obtain ^0Uj/3t^ give 
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Equation (11) and the relations 


b i ^i 



rn m 
2<*i 


a? 111 

1 



+ 



have been used to simplify equations (20) and (21). These relations were obtained by dif- 
ferentiating equation (7) with respect to time t, setting t = 0, and then substituting the 
expressions given by equations (8) and (15) for and ^3u./3t^ , respectively. 

Similarly, 0 0 
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and where 
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The additional relations obtained from differentiating the equation of continuity and em- 
ployed to simplify equations (23) and (24) are 
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Substituting equations (10), (15), (16), and (19) to (24) into equation (8) gives the time 
evolution of as a function of the coordinates in the moving system for small and mod- 
erate times. However, we are more interested in the space-averaged value u^. If the 
product UjUj is averaged (integrated) over space, the integrals of the various terms are 

zero except for those for which the exponential factor reduces to 1, so that the result is, 
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for terms up to order t , 
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where the coefficients of the powers of t are given by 
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Although u.m has been obtained by integrating u.Uj with respect to coordinates in the 
moving system, its value is not dependent on carrying out the integration in that system 
and will be the same in the stationary system. 


RESULTS AND DISCUSSION 

As in the case of nonlinear decay without a mean velocity gradient (ref. 3), the non- 
linear terms in the Navier -Stokes equations here produce a proliferation of a new eddies 
or harmonic components at various wave numbers. These effects are contained in the 
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double, triple, and quadruple summations in equations (15), (19), and (22). In the pres- 
ent case the presence of terms containing mean gradients produces harmonic components 
in addition to those in reference 3. For instance, in equation (17), the nonlinear term 
u i ( 9u k/ 9t ) Q contains mean gradient effects through equations (15) and (16). 

Figure 1 compares the results of the present analysis for dUj/dxg = 0 with those 
from reference 3 for the same initial conditions. Note that these initial conditions sat- 
isfy continuity and give nonzero interaction terms, even for S = dU , / dx 0 = 0 
/ 110? 12 *> i \ vi 

(a^ q? = a^qr 1 = 0, a^q? * 0, and ar'q^ * 0). For S = 0, the only essential difference 

between the two analyses is that the present analysis utilizes a Taylor series, whereas 

that in reference 3 results in a series of exponentials. The first, second, and third 

approximations for the Taylor series correspond, respectively, to retaining terms in 
2 3 

t, t , and t in equation (25). Comparison of the various approximations for the Tay- 
lor series solution and of the results for the Taylor series with those for the exponen- 
tial method should give an idea of the accuracy of the results. 

The results in figure 1 for the present Taylor series solution appear to be accurate 
to a dimensionless time of about 0. 003, whereas those for the exponential method are 
good for much larger times. The Taylor series should still give a good indication of 
the effect of shear on a disturbance, particularly for large shear, where large effects 
of shear at small times might be expected. 

For one accustomed to thinking in terms of turbulence, it may seem surprising that 
the shear component UjUg is not zero for a zero mean velocity gradient. It should be 
recalled, however, that the present disturbance is for a limited number of wave-number 
vectors. Even for turbulence, it is necessary to integrate over all directions in wave- 
number space to obtain zero shear stress for zero velocity gradient. 

The effect of shear on the kinetic energy UjU./2 is shown in figure 2, where u.u^/2 
is plotted against dimensionless time for a low, an intermediate, and a high velocity 
gradient. (The same initial conditions are used for figs. 1 to 4. ) For the low velocity 
gradient S = 20, the energy decays, as for zero velocity gradient. For the intermediate 
velocity gradient, the energy decays at a slow rate at small times and then appears to 
increase. For the large velocity gradient the energy in figure 2 increases monotonically. 

Figure 3 shows the behavior of the individual components of the disturbance for a 

o 

high velocity gradient. All the components except Ug increase with time. Inasmuch as 

production terms (terms multiplied by 6-j) do not occur in the equations for u|, one 

might suppose that component increases because of energy transferred into it from u^ 

by pressure terms (terms with qP in the denominator). Figure 4, where the pressure 

* 2 

terms have been neglected, indicates that is the case, since u| without the pressure 
terms decreases with time. The shear component still increases, apparently because 
the u^ in -UjU 2 increases sufficiently rapidly to offset the effect of the decrease 
in u 2 . 
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All of the results thus far were for the same initial conditions. In order to get an 
idea of how changes in initial conditions can affect the components of the disturbance, 
the directions of the intensity and wave-number vectors relative to the direction of the 
velocity gradient are varied in figures 5(a) to (e). Those figures, together with fig- 
ure 3, present the results for a set of six permutations (out of a possible 36) of the com- 
ponents of the intensity and wave-number vectors for fixed magnitudes of the vectors. 
Comparison of the various figures shows that changes in the initial conditions affect the 
evolution of the disturbance both quantitatively and qualitatively. 

Perhaps the_most interesting aspect of the curves in figure 5 is that, although the 
equations for u^ contain production terms, u^ for dUj/dXg * 0 decays faster than it 
does for dUj/dXg = 0 in several cases. That result might be caused by a change in 
sign of the production term, in which case the shear component u-jUg would become 
positive, or by large energy transfer out of u^ into the other components by the pres- 
sure terms. The results show that both of these effects can be important. In figures 
5(a) and (d) u-jUg does change sign and become positive, so that the production term 
becomes negative. In that case u^Ug tends to increase S rather than to decrease it, 
as it normally does for turbulence. That is, it acts like a pump rather than a brake on 

O 

the fluid. There is also an effect of pressure terms in transferring energy into Ug 
from other components, since u| increases with time even though the equations for that 
component do not contain a production term. 

o 

In figures 5(b) and (c), on the other hand, the abnormally high decay rate of u| is 

due entirely to pressure forces since UjUg is negative and thus the production term 

for u? is positive in those figures. The large effect of directional energy transfer in 
A o o 

those figures is also shown by the fact- that both uJjj and Ug increase with time, even 

though the equations for those components do not contain production terms. 


CONCLUDING REMARKS 

The nonlinear terms in the Navier-Stokes equations produce a proliferation of new 
harmonic disturbances at various wave numbers. The presence of a mean gradient pro- 
duces harmonic components in addition to those produced when it is absent. Although 
the present Taylor series solution seems to be limited to shorter times than is an anal- 
ysis in which exponentials are obtained, it is capable of showing the effects of a mean 
velocity gradient on the early evolution of a disturbance. As the mean gradient in- 
creases, the rate of decay of the kinetic energy of the disturbance decreases. For large 
gradients the energy can increase with time. However, for the chosen initial conditions, 
at least one of the directional components always decayed. 

Rotating the intensity and wave-number vectors of the initial disturbance changed 
the evolution pattern of the disturbance both quantitatively and qualitatively. For some 
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orientations the shear component of the disturbance changed sign, so that it tended to 
increase the velocity gradient (pump the fluid), rather than to decrease it (brake the 
fluid). Also, the shear components of the disturbance could be nonzero even when the 
velocity gradient was zero. These results differ from those for homogeneous turbulence 
apparently because the disturbance considered in this report consisted of a limited num- 
ber of harmonic components. As in turbulence, the pressure forces played a significant 
role in the directional distribution of the disturbance energy. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 13, 1973, 

502-04. 
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APPENDIX -SYMBOLS 


a? 1 defined by eq. (8) or ( 10) 

a? 111 defined by eq. (16b) 

a mnr defined by eq. (20c) 

a mnrs defined by eq. (23d) 

bP defined by eq. (16a) 

b? 111 defined by eq. (20b) 

binnr defined by eq. (23c) 

cf 1 defined by eq. (20a) 

c™ n defined by eq. (23b) 

d? 1 defined by eq. (23a) 

p dimensionless pressure, ( y^jpv^i * 

p* pressure 

q? 1 component of initial dimensionless wave number vector Xgq* 11 

q^ component of initial wave number vector 

q 111 initial dimensionless wave number vector Xgq™* 

q” 1 * initial wave number vector 

S dimensionless mean-velocity gradient, dUj/dxg 

t dimensionless time, pt*f Xq 

t* time 


U i component of mean velocity 

u i component of dimensionless spatially fluctuating velocity 

Uj component of dimensionless velocity, XqU*/i> 

u* component of velocity 

Xj dimensionless position coordinate relative to moving observer, x . - SjjU^XgH 

Xj dimensionless position coordinate, x?/xq 
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Xq characteristic length 

a m defined by eq. (21a) 

o- mn defined by eq. (21b) 

« mnr defined by eq. (21c) 

defined by eq. (24a) 

/3 mn defined by eq. (24b) 

/3 mnr defined by eq. (24c) 

j 3 mnrs defined by eq. (24d) 

6^ Kronecker delta, 1 for i = j, 0 for i * j 

v kinematic viscosity 

7r m definedbyeq. (14a) 

n mn defined by eq. (14b) 

p density 

Superscripts: 


dimensional quantity 

averaged velocity or stationary coordinate 
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Figure 1. - Comparison of present results for a Taylor series with those from ref- 
erence 3 for an exponential method for dimensionless velocity gradient of 0. 
a} - 30(-l,l, -2); af = 120(-2,1, -1); q! = (2,4,1); q? = 2(1,4,21. 
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Figure 3. - Evolution of disturbance components for dimen- sionless velocity gradient of 1000 with pressure terms 

sionless velocity gradient of 1000. aj = 301-1,1, -2); omitted. a| * 30H.1, -2); a? = 1201-2,1, -1); qj - (2,< 

a? ■ 120 (-2,1, -1); q? • (2,4,1); q? = 2(1,4, 2>. q? = 2(1,4, 2). 
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